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Abstract 



In this work, we update the constraints on tree level FCNC cou- 
5^ I plings in the framework of a theory with n isosinglet vector-like down 

quarks. In this context, we emphasize the sensitivity of the B — > 
J/ipKs CP asymmetry to the presence of new vector-like down quarks. 
This CP asymmetry, together with the rare decays B — > Xg ,ill and 
K — > t:vi> are the best options to further constrain the FCNC tree 
level couplings or even to point out, in the near future, the possible 
presence of vector-like quarks in the low energy spectrum, as suggested 
by GUT theories or models of large extra dimensions at the TeV scale. 



1 Introduction 



The study of flavor changing neutral currents (FCNC) in particle physics 
phenomenology has played a key role in the advance of high energy physics 
in the past decades. Already in the early seventies the non-observation of 
FCNC was used to predict the existence of the charm quark through the GIM 
mechanism well before its direct experimental discovery. Subsequently, 
a precise prediction of its mass was made from FCNC kaon processes [Q]. 
Later on, the discovery of the bottom quark and the measure of the B B 
mass difference indicated the presence of a heavy top quark. Due to this 
GIM mechanism, FCNC in the SM arise only at higher loop level and are 
suppressed by powers of light quark masses and small mixing angles. This 
strong suppression make FCNC phenomena a privileged ground to search for 
signs of new physics beyond the SM. In many extensions of the SM there 
is no analog of the GIM mechanism to protect FCNC processes and hence 
potentially large new physics contributions can be expected. For instance, 
a minimal supersymmetric extension of the SM with generic soft breaking 
terms gives rise to dangerously large contributions to FCNC and CP violating 
observables through loop contributions with the supersymmetric particles 
running in the loop . Even more challenging (or perhaps more dangerous) 
than these loop induced FCNCs is the inclusion of tree level FCNC couplings. 

Indeed, the presence of tree level FCNC is phenomenologically, a very in- 
teresting possibility. A minimal extension of the SM with the only addition of 
an extra isosinglet down quark in a vector like representation of the SM gauge 
group induces FCNC couplings in the Z and neutral Higgs boson couplings. 
Some, more realistic models, include at least one vector-like quark per gen- 
eration. These models naturally arise, for instance, as the low-energy limit 
of an Eq grand unified theory ||^. Moreover, vector-like quarks can appear in 
models of large extra dimensions. The existence of extra spatial dimensions 
of the TeV size implies the existence of towers of Kaluza Klein modes for 
the particles that propagate in the new (extra) dimensions. Therefore, if 
some SM fermion propagates in the extra dimensions, for each of these chiral 
quarks and leptons there is a tower of vector-like fermions whose separation 
is of order 1/R, being R the radius of the extra dimensions 0, ^. From a 
more phenomenological point of view, models with isosinglet quarks provide 
the simplest self-consistent framework to study deviations of 3 x 3 unitarity 
of the CKM matrix as well as flavor changing neutral currents at tree level. 
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Of course, these FCNC couplings are severely restricted by the low energy 
results on the different FCNC processes. Nevertheless, it is well known that 
even fulfilling these strong constraints these couplings can have large effects 
on B factory experiments on CP violation [§, |10|, |ll[]. In a recent paper [|ll 



we showed that the possible mismatch between the SM expectations and the 
measured value for the CP asymmetry in the B — > J/ipKg decay can be 
easily explained in a model with an additional vector like down quark. In 
this paper, we complete the previous analysis and generalize it to the case 
where a tower of vector-like down quarks (VLdQ) is present. In first place, we 
update the strong low-energy constraints on the tree- level FCNC couplings 
and then we concentrate on the CP violation observables in the B factories 
as the most efficient observables to find or at least constrain these tree level 
FCNC couplings. Finally, we compare the simplest case with a single VLdQ 
with models with several vector-like down quarks. 



2 FCNC in the presence of isosinglet quarks 

As explained in the introduction, vector-like fermions appear in different 
extensions of the SM, as for instance GUT models or, remarkably enough, 
theories with large extra dimensions on the TeV scale. All these models have 
several new vector like fermions that mix with the SM fermions and can have 
important effects on the low energy FCNC phenomenology. In this paper, 
we will study FCNC processes with external down quarks, and so we will be 
mainly interested in models with isosinglet vector-like down quarks (VLdQ). 
Note that the presence of extra vector-like down (up) quarks generates FCNC 
in the down (up) couplings to the Z and H. In view of these new possibilities 
of having VLdQs in the low energy spectrum, and especially, to reach a rough 
understanding of FCNC effects in models with large extra dimensions ^], 
we present our bounds in a general framework with n additional VLdQs. A 
general analysis with both up and down isosinglets and isodoublets will be 
presented elsewhere 

From the low energy point of view, the model we analyze here has been 
thoroughly described in Ref. ^. Its main feature is that the additional 
VLdQs mix with the three ordinary quarks and consequently the mass matrix 
in the down sector is now (3 + ra) x (3 + ra). This matrix is diagonalized by 
two (3 + n) X (3 + n) mixing matrices and only the left-handed rotation is 
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observable giving rise to the CKM matrix. As we have aheady announced 
in the introduction, the fingerprint of this model is that it allows tree level 
FCNC in the Z and H vertices. To see how these FCNC couplings appear, 
we can work in the basis where the up quark mass matrix is diagonal. Here 
the down quark mass matrix is diagonalized by a (3 + n) x (3 + n) unitary 
matrix, Vap. In this model charged current couplings are unchanged except 
that VcKM is now the 3 x (3 + n) upper submatrix of V , and at low energies a 
non-unitary 3x3 mixing matrix appears at tree level. However, the mixing of 
doublet and singlet weak eigenstates modifies the neutral current couplings 
that in terms of fiavor eigenstates will be V'^ ■ Diag.(l, 1, 1, 0) ■ V 7^ I4X4. So, 
neutral current interactions are given by. 



r — 9 

2^os^ 



ULi Y ULi - dLa Uaf) Y dL(3 
2 siB^Ow J^rr. 



H 



2M, 



w 



ULi m'^ULi - dLa Uai3 mff dLp 
Uap = Ell V:^V,p, (1) 



where Uds, Ubs or Ubd 7^ would signal new physics and the presence of 
FCNC at tree level. 

In the following, we analyze the effects of these new couplings in FCNC 
processes and we update the phenomenological bounds on them. More specif- 
ically, in processes changing fiavor in one unit, AF = 1, that in the SM go 
through electroweak penguin diagrams of order GFaVf*Vtj, we simply con- 
sider the dominant tree level FCNC contributions, order GpUij- Similarly 
in AF = 2 processes mediated by boxes in the SM of order Gpa {V^jVtjY 
we include two different additional contributions, a pure tree level diagram 
with two FCNC vertices, order GpU^j and a mixed SM vertex contribution 
with a tree level vertex, roughly order Gpo: VuVij Uy [|r^. Other new phy- 
sics contributions in this framework will always be suppressed by additional 
loop factors or higher powers of the FCNC couplings with respect to the 



contributions considered here [ 13 1. Therefore the effective low energy theory 
we consider is then identical to the SM with a non-unitary CKM matrix and 
the presence of tree level FCNC as shown in Eq. (|I|). 
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3 Experimental constraints on the extended 
mixing matrix 

As we have seen in the previous section, all the new physics effects in our 
model are encoded in the extended mixing matrix that gives rise to tree level 
FCNC in the Z and H vertices. The minimal extension from the SM would 
consist in the addition of a single vector-like down quark. In this situation 



the mixing matrix can be parametrized as |T4 



V = Ru (Ou, 0) i?24 {O24, 03) Ru (^14, 02) ^f^M (^12, ^13, ^23, 0l) , (2) 

where V^^^j (6*12, ^^13, 6*23, 0i) is 4 x 4 block diagonal matrix composed of 
the standard CKM [|15|, and a 1 x 1 identity in the (4,4) element, and 
Rij {6ij, 0fc) is a complex rotation between the i and j "families". Note that, 
in the limit of small new angles, we follow the usual phase conventions. 
In fact, this 4x4 mixing matrix can represent a good approximation to 
more complete models with several vector-like generations if the mixings are 
hierarchical, similarly to the standard CKM matrix. 

On the other hand, given that the deviations from unitarity of the CKM 
mixing matrix are experimentally known to be small, it is possible to use 
an approximate parametrization of the mixing matrix valid for an arbitrary 
number of vector-like quarks |^ . This is an extension of the usual Wolfenstein 
parametrization of the CKM [T^ in terms of. 
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with {Usd,Ubd,Ubs) general complex numbers and D| = ^27=4 l^iP positive 
real numbers. It is possible to obtain the remaining elements of the 3x3 
submatrix corresponding to the SM mixing matrix as a function of (A, A, /j,, 
0, Dj, Dl, Dl, Usd, Ubd, Ubs)- In fact, taking into account that [Dj, Dl, Dl) 
can be at most of order A^ and (^7^^, Ubd, Ubs) are experimentally constrained 
to be < 0{X^), we can keep quadratic or linear terms in {D^, Dl, Dl) and 
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{Usdi Uhdi Uhs) respectively and we obtain to 0{\ 

i2 \4 1 , c /l2,,2 n2 n2 on2 



ud 



^ y \- 1 + 8AV Dj Di-2Di 
2 8 16 2 4 

^ + \Re{Usd} + 0{X') 

o 

,r _ , A' 1 + 4A%4 , 1-4A' + 16AVcos0,6 

V^ro — i — — — A -p A — 

2 8 16 

_ f^A^ + AA2Re{f/,J - ^ + O(A^) 
z 4 o 

2 2 2 8 ^^ 

/ 42 \ n2 _ n2 

= -A + f ^ - A V-^*j + ^^V^^ + ^^"^ + 
- A^ + 0{X') 

Vtd = A(l-/ie-*^)A3 + ^/xe-^<^A5 + ^D2^3 + f/5rf + 
I {Dl - Dj) (1 - ^^e-^^) X' - AX'Usd + O(A^) 
Vts = -AX'+(^-Ai^e-^Ax' + U,, + ^(D',-Dl)x' + 



|a6 + 0{X') (4) 

This is a completely general parametrization of the 3x3 submatrix of a 
complete (3 + n) x (3 + n) unitary matrixB. Hence, it includes also the 
simplest case of a single vector-like quark. In this last case, there are several 
relations among the 13 parameters of this general matrix and only the nine 
independent parameters of a unitary 4x4 mixing matrix remain. In fact we 
have, 

Dl = \V,d\', Dl = 1^4.1^ Dl = \V,,\', 
Usd = -VIVm, Uts = -V:,V^s, Um = -V:,Vu. (5) 



^Notice that this parametrization is approximate to 0{\^) but an exact solution can 
be obtained numerically 
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Hence, we have 4 relations among the parameters in Eq. 

\UJ' = DlDl W = DlDl 

\Um? = DIdI UsdJbsU;a = DlDlDl (6) 

In the general case these equalities are replaced by inequalities, 

\Us,?<DlDl \U,s?<DlDl 

\Um? < DlDl (7) 

Naturally, the elements of the extended mixing matrix corresponding 
to CKM elements are experimentally measured in low energy experiments. 
Some of these measurements are obtained from tree-level processes, therefore 
not affected by new physics to a high degree of approximation. Hence they 
constrain directly the corresponding elements of our extended mixing matrix. 
Specifically, we have p|, pi, at 95% C.L., 



• 0.2150 < \Vus\ = A < 0.2242, 

• 0.0364 < iKfel = < 0.0440, 

• 0.074 < iKb/Kfel = /iA < 0.106, 

• 0.9719 < \Vud\ < 0.9751, 

• 0.192 < \Vcd\ < 0.256, 

• 0.948 < \Vcs\ < 1.0. 

These constraints restrict the different parameters in Eq. (H), and so from 
the bounds on \Vud\ and \Vcs\, we obtain < 5.5 x 10~^ and Dl < 5.6 x 10"^ 
respectively. Another constraint |^ H, ^ comes from the SU{2)l coupling 
of the Z° to hb. In the SM, this coupling is (V^t^^j ■ VcKM)bb = 1, but in this 
model it is modified to Uhb] hence, we have [|l9l -D^ < 9 x 10~^. 

Notice that, in the general case, the parameters are completely inde- 
pendent from the FCNC couplings Uap- However, in a more definite model, 
as for instance the single vector-like quark model, due to the unitarity of the 
4x4 matrix, these constraints have a strong impact on all other elements of 
the extended mixing matrix and consequently also on the tree level FCNC 
couplings, as shown in Eq. (Bf). In any case, even with these restrictions. 
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the FCNC couplings can have large effects on rare processes, where the SM 
contributes at the 1 loop level and new physics is allowed to compete on 
equal footing. In fact, most of the experimental results are well accommo- 
dated within a pure 3 generations SM and hence these measurements provide 
additional constraints on the FCNC couplings. In other observables, like 
CP asymmetries or B rare decays, the experimental results may still differ 
from the SM prediction once the experimental accuracy is increased in the 
near future. Therefore, we analyze these observables from a slightly different 
point of view, and in the last section we show what are the possible deviations 
from the SM consistent with the constraints discussed here. 

In the first place we analyze the constraints associated with kaon physics 
and then the constraints we get from the B system. 



3.1 Kaon physics constraints 

Rare kaon decays and CP violating observables in the kaon sector can re- 
ceive important contributions from the FCNC coupling Usd- In fact, this 
new coupling is constrained mainly by the decay [K^ A*^/^ )5d and the 
experimental value of e' /e. Other observables constraining this coupling are 
K+ TT+z/z/ and Ek H- 

The decay Kl — >• /x"*"//" is CP conserving and in addition to its short dis- 
tance part, given by Z penguins and box diagrams, receives important contri- 
butions from the two photon intermediate state which are difficult to calcu- 
late reliably. Unfortunately, the separation of the short-distance part (similar 
to — > vr+i/z/, free of hadronic uncertainties) from the long-distance piece 
in the measured rate is rather difficult. Therefore, the full branching ratio 
is generally written as a sum of a dispersive and absorptive contributions, of 
which the latter can be calculated using the data for Kl — > 77. The disper- 
sive contribution can be decomposed as a long distance and a short distance 
part. Following [^D|, we can write down. 



Br {Kl ^ iHi)sD = 6-32 x IQ-^ [Cu2Z Re {Usd) 

+ A, + Fo(xt) Re(Af)]' < 2.8xl0-^ (8) 

where Cu2Z = -{V^GfM^/ti^)-' ^ -92.7, Af = V*,Vi,, A, = -6.54 x IQ-^ 



is the charm quark contribution and Yq is the Inami-Lim |22] function as 



7 



defined in [EH]. Experimentally, we must require this branching ratio to be 



BR{Kl H+fi-)sD < 2.8 X 10-9 IT 



A second observable constraining Ugd is e'/e. Within the SM there is 
a cancellation between QCD and electroweak penguin contributions (domi- 
nated by Z penguins) which suppresses this ratio. When new physics enters 
into the game and if, as it is expected on general grounds, it does not af- 
fect considerably the QCD contributions but does so with the Z penguins, 
the abovementioned cancellation does not take place and significant devia- 
tions from the SM results (or strong bounds on the new parameters) can be 
expected. We now decompose e'/e as follows. 



(3uCu2zlm{Usd) + Alm(Af) 
(3u = [1.2 - Rs rz 5f ^'^] 



A = /5c7-Co - 2.3 + Rs [1.1 4'^'^ + 

(1.0 + 0.12 rz) i?f '^]. (9) 

where the first term comes from the Z piece and the other one contains all 



the remaining ones ||20|. It is worth noting that unlike previous observables, 
here the theoretical errors overwhelm the experimental precision. The main 
sources of uncertainties lie in the parameters, B^^'^^ and B^^'^\ The impor- 
tance of these uncertainties is somehow increased because of the cancellation 
we have mentioned. 

Here, we take i?, = 1.5 ± 0.5, rz = 7.5 ± 1 and Bf^^ = 0.8 ± 0.2 ^ 



The value of B^^"^^ has caused some controversy in the literature because 
of the different values obtained in different schemes. We take two different 
values in order to illustrate the situation where new physics in the s-d sector 
is needed or not needed to reproduce the experimental value of e'/e. For set 
I, we use B^^'^^ = 1.0 ± 0.2 with the other parameters as given above, as 



in Refs. [pO| , p^ , |25| which comes from large Nc calculations [^ and lattice 
analysis, and this tends to favor the presence of new physics in Uds- In 
set II, we use B^^"^^ = 1.3 ± 0.5 in order to incorporate the predictions of 
Refs. [^, where inclusion of the correction from final-state interactions in 
a chiral perturbation theory analysis tends to favor the SM range. Still, there 
are other schemes where different values for Bq^^^^ and B^^"^^ are obtained 



2^, pn, pTl. In fact, we have explicitly checked that with the values in [p£ 



B^^^^^ = 2.5 ± 0.4 and sF^^ = 1.35 ± 0.20, the constraints on the FCNC 
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couplings Uai3 are still consistent with the results presented below. 
Numerically, we have for the central values 



e 



-7.8 Cu2zlm{U,d) + 



6.2 
9.9 



Im (^Xf) 



(10) 



with the parameters in set I and set II respectively and we calculate its errors 
with a Gaussian method. This observable has to reproduce the experimen- 
tally obtained value of {e'/eY^'P = (2.12 ± 0.46) x 10"^ 

A theoretically cleaner constraint although less restrictive than the 
previous two constraints is provided by BIi{K~^ tt'^uu), 



1.55-10- 



Cu2zIm{Usd} + Xo{xt)-ImXf) (11) 



+ Cu2zRe{U,d} - 2.11 ■ 10-^ + Xo{xt)ReXr 



with Xo(xt) = Co{xt) — ABQ^Xt) a gauge invariant loop function combination 
of boxes and vertices This decay has already been measured in the 



experiment E787 at BNL |Q, however, so far, a single event has been found 
and this is not enough to provide a definite value for the branching ratio. 
Hence we take here only the upper limit at 95% CL. BIi{K^ n^uu) < 
8.3 X 10-1° i. 

Finally, we include also ek, whose leading-order expression is [^, 0: 



Im <^ -7]tt {Usdf + 



a 



in sin 9w 



l,]=C 



(12) 



where 5*0 is another Inami-Lim function and the QCD correction factors 
(which take into account short distance QCD effects) are given as follows. 



1.38 ±0,20, 



rjtt 



0.57 ±0.01 



rjtc 



0.47 ±0.04 



(13) 



Here, contrary to Ref. PDI, the coefficient Iq (2^) of the linear term in Uds is 
characteristic of the present model, therefore, in principle, the irrelevance of 
Sk to constraint Uds is not fully guaranteed. 



^In Ref. the old value of this BR at 1 cr was used, therefore both a lower and an 
upper limit were considered 
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Figure 1: Effect of tlie constraints from Ek, n^uu, Ki fi^jj, and 

e'/e on tlie Ugd FCNC coupling. 



At tliis point, it is important to empliasize tliat Eqs. (|TT]) and 

(|12D are completely valid in the general model with n additional VLdQs. 

In Fig. we present the effects of the different constraints in the Usd 
coupling. As we can see in this figure, the most efficient constraints are 
provided by Kl — >■ that constrains Re(Usd), and e'/e, that bounds 

lm{Usd) to the rectangular box in the center. However, these constraints will 
be difficult to improve due to the large hadronic uncertainties. Similarly, 
the constraints from ex are very precise on the experimental side and the 
precision of this limit on Usd is determined by the hadronic parameter Bx = 



0.85 ± 0.15 |23, |25|. Hence, it is not expected to improve largely in the near 



future. On the other hand, the decay — > Tr^iyiy is much cleaner from 
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the theoretical point of view [Q. The experimental value for the branching 
ratio, based in the single event found so far is, 



BR{K+ TT+uu) = 1.5 



+3.4 
-1.2 



X 10 



,-10 



(14) 



which gives an upper bound of. 



BR{K 



■+ 



vr+z/z/) < 8.3 X 10 



-10 



at 95% CL. 



(15) 



Unfortunately, it is clear from the errors that we cannot obtain a lower bound 
on this BR at 80% CL, still we have included in the figure a dashed circle 
showing the effect of a future lower limit which would correspond to a hy- 
pothetical value of 8 X 10~^^ and would exclude the region to the left of the 
small rectangle. Moreover, it is interesting to notice that improving the up- 
per bound to the level of 2 x 10~^° would provide a bound on the same level 
as the combined bounds from Kl — > yU^yU^ and e'/e. In fact, these values 
could be reached after the analysis of the stored data from E-787 and the 
sensitivity of the new experiment E-949 fS^ (already approved) will reach 
10~^^ in the next few years. Hence, this decay will provide the most stringent 
and clean bounds on the Usd coupling in the near future, although at present 
the bounds are still obtained from — > and e'/e. 

In Figs. D and |], we present a scatter plot of the allowed values of Usd 
in the general model with n additional vector-like down quarks, although we 
must emphasize that this allowed region does not change at all even if we go 
to the more restricted case of a single VLdQ. We impose all the constraints 
described above, with the e'/e parameters from set I and set II in Fig. Q and 
^ respectively. In both approaches the bound for the real part is. 



This constraint is directly obtained from Eq. (P) with the limit values of 
Re(Af ) which are Re{Xf)max = -4.9 x 10"^ and Re(Af )™„ = -2.1 x 10"^ 
from Eq. (^ to 0{X^) for = vr, respectively. Therefore, this implies that 
the additional correlations among the parameters in the 1-VLdQ model are 
irrelevant for the Usd bound. 

The constraints on the imaginary part of this coupling depend slightly 
on the adopted value for Bq'^'^\ For set I, with Bq^^'^^ = 1.0 ± 0.2, we can 
see that Usd is necessarily positive and does not reach the origin, indicating 



1.3 X 10"^ < Re{Usd} < 4.0 x 10 
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(16) 
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0.5 

10^ Re (U 



Figure 2: Phenomenological bounds on the FCNC coupling Ugd with the e' /e 
hadronic parameters in set I. 
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0.5 

10^ Re (U 



Figure 3: Phenomenological bounds on the FCNC coupling Ugd with the e' /e 
hadronic parameters from set II. 



13 



marginally the need of new physics for e' je. The allowed range is 1.0 x 
10"^ < \m.{Uds\ < 6.5 X 10~^. On the other hand, in set II, we get an 
allowed area, —1.5 x 10"^ < lm{Uds} < 6.0 x 10"^, including the SM as 
one of the possible points in agreement with the experimental results. This 
bounds are directly obtained from Eq. (|10D with lm{X^'^)max = 1-4 x 10~^ 
and Im(A^'^)mm = 6 x 10~^ assuming a 50% error, which corresponds to the 
Gaussian errors that we use in this constraint. In view of these two different 
options to calculate e'/e, we take a conservative bound on lm{Usd) including 
both possibilities, 

-1.5 X 10""^ < lm{Usd} < 6.5 x 10~^ (17) 

It is interesting to notice, that these bounds turn out to be much more strin- 
gent than the direct bounds on FCNC couplings usually quoted in the litera- 
ture for VLdQ models [0. This improvement is mainly due to the inclusion 



of the SM contributions that were completely neglected in the calculation 
of Ek in the bounds presented in the improvement of the experimental 
results in {K^ — > fi''^fJ'~)sD, K'^ tt+z/z/, and to the inclusion of the bound 
from e'/e. In this way, our bounds basically agree with the general bounds 
in pg. 



3.2 B physics constraints 



In this section, we study the constraints on the Um and Ubs couplings from 
B physics FCNC. In fact, the experimental information in the B system is 
improving rapidly with the new results from B factories and so, it is impor- 
tant to update the bounds on these FCNC couplings. In all the following 
processes, the choice of set I or set II to calculate e'/e has no relevant effects 
on the Ubd and Ubs couplings. Hence, we analyze the constraints making no 
distinction on the hadronic parameters used to calculate e'/e. In first place, 
we concentrate on the constraints from CP conserving processes and then we 
add the information from the B^ CP asymmetries. The main CP conserving 
processes constraining Ubd and Ubs are BR(i? — > Xd,sl^l~) and AMb^ 

From the upper bound on BR(i? 
suming BR(i? 



Xdl+l-) < BR{B 



Xsl 



Xsl- 
l- 



1-) < 4.2 X 
) we have 



10-^ 36 



and as- 



T[B 



X,V 



v 



0, 

U2Z Ubq 



r 



TT^sin^^ 



w 



\Vcb\ firric/mb) 



X 
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Figure 4: Constraints on the Um coupling from Bd-Bd mixing in the n-VLdQ 
model (left) and in the 1-VLdQ model (right). The circle shows the bound 
obtained from the B — > Xdl^l~ 



2 sin^ e 



w 



sin^e 



w 



-j < 4 X 10"^ (18) 



with f{z) 



Iz'^ + 8^6 



24z^ log z) a phase space factor due to the 



mass of the charm quark. From here we get, 

Fo [xt] + Cu2z Ut, 



Replacing here A^^ 
get a bound \Ubd 



AX^l- fie 



<0.15. (19) 

- 8 X 10-3 (1 - 0.4e-*'^) with 



1, it is clear that the A^ contribution can be safely neglected and we 



< 1.7 X 10" 



However, this is not true in the case of 



A 



bs 



-AX' 



0.040 that shifts the \Uh«\ constraint to a circle of radius 



1.7 X 10-3 centered in (-0.040/92.7,0). 



Additionally, from B B mixing we have |Ty, [TT], |T2 



fit 



71 



127r^ 



S'o(xt) Ai^^Abg, 



A 



hq 



1 



3.3 



Ubq 
X'" 



+ 165 



(20) 
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10^ Re (UJ 



10-^ Re (UJ 



Figure 5: Constraints on the Ubs coupling from B 
model (left) and in the 1-VLdQ model (right). 



xj+i- 



in the n-VLdQ 



23 , and the mass difference 



where the new parameters are defined in Ref. 
AMb = 2 |M]^2' I • The experimental values for these observables are AM^ , 



(0.472 ± 0.017) X 10^2 g-i and AMb, > 14.9 x 10^^ g-i 95% C.L. 

At this point, it is important to compare the bounds that can be ob- 
tained from the mass difference and the semileptonic decay. From BR(5 
Xd,sl^l~)-, we get a bound on |Aj'^Fo + Cu2zUbq\'^ in Eq. (p!9[), while from the 
mass difference, neglecting the linear term, we obtain a constraint on the com- 
bination |(Ai^)^S'o + ^Cu2zUlq\. Hence, it is clear that the relative size of the 
tree level FCNC with respect to the SM contribution is always bigger in the 
semileptonic decay due to the presence of an additional factor Cu2Z — —92.7. 
Nevertheless, the experimental constraints from BR(i? — > Xd,s^^l~) are still 
much larger than the typical SM prediction while the AMb^ measurement 
is already saturated by the SM contribution. This implies that the AMb^ 
constraint is more effective in the case of Uhd and both experiments give rise 
to constraints of the same order of magnitude, dominating at the end the 
AMb^ bound. This can be seen in Fig. ^, where we show the allowed region 
of the parameter space both in the general model with an arbitrary number 
of VLdQs and in the minimal model with a single VLdQ. In these figures. 
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the constraint from BR{B Xdl^l ) is shown as a circle shghtly shifted 
from the origin and a radius of 1.7 x 10~^. However, as we can see here, 
the bounds on Um are |Re{[/ferf}| < 0.7 x 10"^ and \Iui{Um}\ < 1-2 x lO'^, 
which are directly obtained from the AM^^ constraint as the envelope of the 
curves liXffSo + 4Cu2zUid\ with all the allowed values of . Notice that 
already in this most constrained model, i.e. the model with a single VLdQ, 
the Uhd coupling is only bounded by these processes and therefore we obtain 
the same bounds both in this case and in the general case with n VLdQ. 

In the case of the Ubs coupling the situation is completely different. In 
first place AMb^ has not been measured yet, and so only a lower bound on 
the mass difference is available which is not useful to set a constraint on Ubs- 
Moreover, it is easy to see from Eq. ( pUD that for similar values of Um and Ubs, 
the FCNC effects on AMb, will be suppressed by a factor (Af'^/Af'^)^ ~ A^ 
when compared with the effects on AMb^- This implies that in this model 
we cannot expect to observe the FCNC effects in Bg-Bg mixing I. Hence, 
this coupling is only constrained by the upper bound on BR(5 — »• Xsl~^l~). 
In Fig. ^, we show the constraint from BR{B —>■ Xsl~^l~) as a circle slightly 
shifted from the origin, that implies an upper bound of \Ubs\ < 2 x 10~^ 



[37, 41 1 both in the general model with n VLdQs and in the model with a 
single VLdQ. Similarly, the b ^ s'j branching ratio provides constraints of 
the same order of magnitude ^2 . 



Despite these strong constraints from CP conserving processes, the CP 
asymmetries in B decays are still very effective to constrain the Ubd coupling 



J3, S]. Recently, the arrival of the first measurements of the B J/ipKs 
CP asymmetry, aj/^, from the B factories has caused a great excitement in 
the high energy physics community. 



These values correspond to a world average of aj/^ = 0.51 ± 0.18, that can 
be compared with the SM expectations of 0.59 < ajj^ = sin(2/3) < 0.82 
with P = aig{—VcdV*ii/(ytdV-i-l)). The errors are still too large to draw any 
firm conclusion. In fact, if we take the world average at 95% C.L. we do 



0.34 ±0.20 ±0.05 (BaBar 
0.58l°:ilH^ (Belle g§) (21) 

0.79l[!1^ (CDF 



•^This is not always true in the presence of an up vector-like quark, as can be seen in 
Ref. ioH 
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Figure 6: Allowed region for the Um coupling requiring a value of aj/^ asym- 
metry to reproduce the world average at 1 o" in the ra-VLdQ model (left) and 
1-VLdQ model (right). 

not get any improvement over the constraints from CP conserving processes. 
However, anticipating the improvement of the experimental errors from B 
factories, we take the world average at 1 cr level to see the effects on the Ubd 
coupling. From Eq. (|20|) , it is straightforward to obtain aj/^ as, 

aj/v, = sin (2/5- arg( Am)). (22) 

Using the world average at 1 a as the experimentally allowed range, we show 
in Fig. 1^ the resulting region for the Ubd coupling. As we can see here, 
the Ubd allowed region is sizeably modified from this constraint. The outer 
regions in the second and fourth quadrants are reduced while the central 
region corresponding to the SM remains filled; this situation represents an 
improvement over the analysis presented in Ref. |^^. Note that the results 
are essentially similar for the n-VLdQ and the 1-VLdQ cases. As we can see 
here, with the world average for the aj/^ CP asymmetry there is no need of 
new physics in the Ubd sector. 
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4 Discovery potential in B factories 



As we have explained above, the present results in the aj/^ asymmetry in 
Eq. (|2l|) are not precise enough to make a definitive statement about the 
presence/absence of large new physics contributions in B decays. Still, these 
measurements, and especially the BaBar value which is the most precise one, 
leave room for an asymmetry considerably smaller than the standard model 
expectations. The SM range is certainly outside the la BaBar range but not 
outside the world average. This potential discrepancy is at the origin of sev- 
eral papers ^] studying the implications of a small aj/^ in the search of 
new physics. The papers in Ref. provide general parametrizations of new 
physics contributions to B mixing and decays. They essentially show that 
this possible mismatch among the measured value and the SM expectations 
may be due to the presence of new physics either directly in B physics or 
in K physics modifying indirectly the unitarity triangle fit. However, these 
papers do not provide a definite new physics example fulfilling this task. On 
the other hand, the two papers in Ref. analyze SUSY models, where 
no sizeable effects in B mixing and decays are generally expected and the 
unitarity triangle fit can be modified only through new SUSY contributions 
to K physics. In the following, we show that a model with tree level FCNC 
from the mixing with vector-like quarks would be a natural candidate for a 
model with clear deviations from the SM expectations in CP asymmetries, 
satisfying simultaneously all other experimental constraints. 

From the expression of aj/^ in Eqs. (|20|) and (|22|) , with lA^'^l G [5 x 
10~^, 1.2 X 10~^] and \Ubd\ ^ 10"^, it is clear that arg(Aferf) can be very large. 
Indeed, in the general model with an arbitrary number of VLdQ, there is 
no further infiuence on Um from constraints on other elements of the mixing 
matrix. In fact, only the inequalities in Eq. can restrict this coupling, 
however, given the bounds on Df, they are almost inoperative in this case. 
Hence, all the points inside the allowed contour in Fig.^ are equally probable 
and any value of the asymmetry is possible with a comparable probability. 
Still, the especial geometry of this allowed region and the preferred orien- 
tation of A^^ slightly favors the SM range over other values. Furthermore, 
even in a more constrained model, as in the model with a single VLdQ large 
departures from the SM range are possible. The expected values for the 
asymmetry both in the n VLdQs and in the 1 VLdQ model are shown in 
Fig. ^ as an histogram of the distribution of 95000 events satisfying all the 
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Figure 7: Distribution of the values of ajj^ for 95000 events in tiie n-VLdQ 
(left) and in the 1-VLdQ (right). 
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Figure 8: Correlation between the values of aj/^ and \Um\ in the model with 
a single VLdQ. 

constraints in the previous sections. As expected from the above discussion, 
in the general case all possible values of the asymmetry have similar proba- 
bility. On the other hand, in the model with a single VLdQ, this distribution 
is clearly peaked on the SM range. Moderate deviations are rather frequent 
and larger departures are more rare but still possible for any value of the 
asymmetry. 

In Figs. 1^ and |^, we show, in the 1 VLdQ model the correlation of the 
possible values of the asymmetry with \Uhd\ and \Ubs\ respectively. Here, 
we see that the range 0.55 < aj/^ < 0.85 corresponding to the SM expected 
range concentrates most of the events and can be reproduced with \Uhd\ = 
and any allowed value of \Ubs\- However, for different values of the asymmetry 
there is a clear correlation between aj/^ and the minimum value of \Um\ 
required to obtain this asymmetry. For instance to obtain an asymmetry 
below 0.5, a \Ubd\ > 2 x 10~^ is needed. This required minimum is also 
true in the general model with an arbitrary number of VLdQs. Similarly, 
we see in Fig. ^ that these large values of Ubd correspond to low values of 
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Figure 9: Correlation between the values of aj/^ and \Ubs\ in the model with 
a single VLdQ. 
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Figure 10: Allowed region for the Um coupling consistent with the BaBar 
result for the aj/^ asymmetry at 1 o" in the n-VLdQ (left) and in the 1-VLdQ 
(right). 



Ubs fTTI) however it is important to emphasize that this correlation is only 
true in the minimal model, but not in a model with an arbitrary number of 
VLdQs. To see this we show in Figs. and |TT] both in the n-VLdQ and in 
the 1-VLdQ models an scatter plot of Re{Ubd) versus lm{Ubd) and Re{Ubs) 
versus lm{Ubs) for an asymmetry corresponding to the BaBar result at la, 



0.14 < aj/^ < 0.54 In this Fig. ^ we see that the great majority of 



the allowed points are in the range 1 x 10~^ (2 x 10~^) < \Ubd\ ^ 1-2 x 10~^, 
in the n-VLdQ (1-VLdQ) model, i.e. a large, non-vanishing Ubd coupling is 
required to reproduce the BaBar asymmetry. In particular, this means that, 
within this model, a low CP asymmetry implies the presence of new physics 
in b-d transitions. This conclusion would be unchanged in the general model. 
On the other hand, we see that, for these points, in the 1-VLdQ model, the 
coupling Ubs is always restricted to the range \Ubs\ ^ 3 x 10^"^; hence all the 
allowed points have simultaneously high \Ubd\ and low \Ubs\- Indeed, it is easy 
to obtain in the 1-VLdQ model, from Eq. (H), the relation UbdUbs = —UsdDl- 
The region in the Uds plane does not change with the inclusion of the aj/^ 
constraint, and then we still have, \Usd\ ^ 6 x 10~^ and Dl < 0.009. Taking 
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Figure 11: Allowed region for the Ubs coupling consistent with the BaBar 
result for the aj/^ asymmetry at 1 o" in the ra-VLdQ (left) and in the 1-VLdQ 
(right). 

into account that a low aj/^ requires \Ubd\ > 2 x 10~^, this clearly implies an 
absolute upper bound, \Ubs\ < 3 x 10^"^. Therefore, for this set of points, we 
can not expect a new-physics contribution in the b ^ s transition. However, 
this relation is not valid in the n-VLdQ model and therefore this correlation 
is lost as can be seen in Fig.|TT|. 

Also in Fig. |10|, we find a few points (~ 0.1% of the points) which have 
simultaneously \Ubs\ ^ 1 x 10~^ and \Ubd\ ^ 3 x 10~^. This second class of 
points is only possible in the vicinity of the SM and they disappear if the 
value of the asymmetry is reduced to aj/^ < 0.52 i. 

Therefore, we can conclude that in a model with a single VLdQ a low 
CP asymmetry aj/^ < 0.52 implies \Ubd\ > 2 x 10~^ and simultaneously 
\Ubs\ < 3 X 10""^. On the contrary, in a general model with an arbitrary 
number of VLdQs, we still have \Ubd\ > 2 x 10~^, but this has no influence 
on the Ubs coupling, and it is only constrained by the B — > Xgll branching 

^ Still, it is important to emphasize that these points also require the presence of new 
physics in B decays. 
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ratio. 

At this point, it is very interesting to examine the predicted branching 
ratios of the decays B X^^s^l for this set of points. From Fig. [1^, where 
we have included the circle corresponding to the experimental bounds in 
these decays, it is clear that we can also expect a very large contribution 
to i? — >■ Xiill. In this case, the branching ratios for the decays are 
strongly enhanced from the SM prediction, reaching values of 1.0 x 10~^ < 
BR{B XJ+l-) < 1.8 X 10"^ and 6.0 x 10"^ < BR{B ^ Xduu) < 1.0 x 
10~^. While, on the other hand, the low values of Ubs imply that the Xs 
decays remain roughly at the SM value. Conversely, in the points with small 
Ubd and large Ubs, there is no sizeable departure from the SM expectations in 
B — >■ XdU and the Xg decays are now close to the experimental upper range. 
Namely, we obtain, for the points to the right of Fig. |lT|, with Re [Ubs) — 1 x 
10~^, BR [B —>■ Xgl'^l^) — 2.7 X 10^^ to be compared with the experimental 
upper bounds of BR{B — * Xgl^l^) < 4.2 x 10~^. However, this possibility 
is only marginal in the la BaBar range for the minimal model with a single 
VLdQ. In any case in the general model, the possibility of large Ubs couplings 
is again open. For analysis of the phenomenological effects of this coupling in 
radiative b s'y decays, Bs-Bs mixing and rare decays in a slightly different 



context, see Refs. m, 42 



5 Conclusions 

In this work, we have updated the constraints on tree level FCNC couplings 
in the framework of a theory with n isosinglet vector-like down quarks. The 
inclusion of the constraints from {K^ ^~)sd-, ^' I^-, and K -nvv has 

allowed to improve sizeably the bound on the Usd coupling. The precise range 
of allowed values for this coupling depends strongly on the hadronic input in 
the calculation of e' je. Our summary in Table |I] includes the main theoretical 
approaches [|0|, ||, ||, l^, |2|, |2|, m. In the near future, 



vr+z/z/ will be quite useful to further constrain Usd without large theoretical 
uncertainties. 

We have calculated the constraints on Um and U^s from AM^^ , B — > X^ Jil 
and the B — >■ J /ipRs CP asymmetry. In Table ^ we summarize our results 
on these bounds without the additional constraint from the B J/ipKs CP. 
Note that all these bounds on Ubq are independent of the hadronic inputs in 
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Re {Usd} 


> -1.3 X 10-5 


< 4.0 X 10-6 


Im {U,a} 


> -1.5 X 10-6 


< 6.5 X 10-6 


\Re{UM}\ 


< 7 X 10-4 


\lm{UM}\ 


< 1.2 X 10"^ 


\Ubs\ 


< 2 X 10"^ 



Table 1: Constraints on the tree level FCNC couplings from rare processes 
in the K and B systems. 



e'/s. In addition, we have shown that the B — > J/ipKs CP asymmetry is 
especially sensitive to the presence of a Ubd FCNC coupling. In this regard, 
we have shown that the 1 a value of the world average is already able to 
exclude half of the allowed region for the Uij,i coupling. However, due to 
the still large experimental errors this constraint has a small effect at 95% 
C.L. Assuming a low value of the aj/^ asymmetry below ~ 0.52 (following 
BaBar range at 1 a) we have shown that, in models with n VLdQs, tree-level 
FCNC in the b-d sector are mandatory. The FCNC coupling Um is required 
to be greater than 1 x 10"^, rising to 2 x 10"^ in the 1-VLdQ model. The 
Uhs coupling in the general model is bounded by 2 x 10"^, although for the 
1-VLdQ model the bound goes down to 3 x 10"^. 

Therefore, a clear favorable scenario for these models would be to find 
simultaneously a low aj/^ (below ~ 0.5) and a large BR(Bd — > XJl), at least 
one order of magnitude bigger than the SM expectations. Values of BR(i?d 
Xgll) close to the present experimental bounds are still possible in the n- 
VLdQ model. Nevertheless, in the 1-VLdQ model these BR are expected to 
be similar to the SM values. In summary, this aj/^ CP asymmetry, together 
with the rare decays B — > Xg^dl and K — > tti/P are the best options to 
further constraint the FCNC tree level couplings or to discover the presence of 
vector-like quarks in the low energy spectrum as suggested by GUT theories 
or models of large extra dimensions at the TeV scale. 
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